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I.  Introduction  ; 

In  [1]  an  adaptive  constant  false  alarm  rate  (CFAR)  detection  algorithm  is 
developed  from  the  generalized  maximum  likelihood  ratio  (MLR)  test  by  using  the 
experimentally  obtained  result  in  [3],  that  most  optical  images  can  be  modeled  as  a 
whitened  Gaussian  random  process  with  a  rapidly  space-varying  mean  and  a  more 
slowly  varying  covariance.  Such  a  CFAR  test  is  closely  related  to  the  test  developed 
by  E.  J.  Kelley  in  [4]  for  detecting  radar  targets.  The  probability  of  false  alarm 
(PFA)  of  the  CFAR  test  in  111  is  a  function  only  of  N,  the  number  of  samples,  and 
K,  the  number  of  reference  image  scenes  used.  The  CFAR  detection  algorithm  in  [1] 
allows  one  to  find  a  detection  threshold  which  achieves  a  fixed  PFA  over  the  entire 
set  of  image  scenes  which  is  invariant  to  intensity  changes  in  the  noise  background. 

The  CFAR  detection  algorithm  considered  in  [1]  is  suitable  only  for  detecting  a 
target  pattern  in  one  main  image  scene  and  a  number  of  other  noise-only  reference 
image  scenes  which  contain  negligible  signal  energy.  However,  in  many  applications, 
one  needs  to  test  for  the  presence  of  a  signal  pattern  which  has  nonnegligible  unk¬ 
nown  relative  intensities  in  several  optical  bands.  As  a  consequence  it  is  of  impor¬ 
tance  to  generalize  the  previous  CFAR  detection  algorithm  [1]  to  a  test  which  is  able 
to  detect  the  presence  of  an  optical  signal  pattern  with  non-zero  intensity  in  several 
signal-plus-noise  bands  or  channels.  A  effort  was  made  in  [2]  to  find  and  compute 
the  statistics  of  this  generalized  CFAR  test,  but  the  results  were  incomplete. 

In  this  paper  the  approach  first  >nsjdered  in  [2]  to  find  this  new  CFAR  test  is 

improved  and  solved.  In  Section  II  tiiv.  general  hypothesis  test  is  formulated  and 

found  in  terms  of  the  generalized  maximum  likelihood  ratio  principle.  The  result  is  a  _ _ 

r 

CFAR  test  for  a  signal  with  unknown  relative  intensities  in  J  channels.  If  J  =  1,  the  — — 
resulting  test  reduces  to  the  standard  normalized  matched  filter  test  for  finding  a  sig¬ 
nal  in  clutter  of  unknown  and  varying  intensity.  The  detection  statistic  found  for  this  ! - 

new  test  is  similar  to  the  adaptive  array  test  for  spread  spectrum  communications _ 

t 

obtained  by  Brennan  and  Reed  in  [5]  except  that  their  test  was  not  obtained  from  a  -  — — 

<■*  y  Codas 

.Avail  and/or 
oist  Special 

!  A  A  I 


□  LI 
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hypothesis  test.  The  test  in  [5]  for  automatic  synchronization  was  derived  as  a  least 
mean  square  criterion. 

In  order  to  analyze  the  performance  of  this  new  CFAR  test,  the  probability  den¬ 
sity  function  of  the  MLR  test  for  both  hypotheses  is  found  exactly  in  Section  III. 
These  probability  densities  are  used  then  to  calculate  the  probability  of  a  false  alarm 
(PFA)  and  the  probability  of  detection  (PD)  as  a  function  of  the  detection  threshold 
in  a  manner  similar  to  that  utilized  in  Til.  The  PFA  of  the  test  is  computed  in  a 
closed  formula  which  is  independent  of  the  covariance  matrix  of  the  actual  residual 
clutter  noise  encountered.  The  probability  density  obtained  here  for  a  real  optical  sig¬ 
nal  in  residual  clutter  noise  is  similar  to  that  found  in  [5]  for  complex  communication 
channels.  However,  the  method  of  derivation  is  different  and  can  be  extended  to 
other  more  complex  detection  problems. 


II.  Formulation  Of  The  Problem 

The  present  detection  problem  is  formulated  in  a  manner  to  that  used  in  [1]. 

First  let  the  column  vector,  {  z(n)  =  [  xx(,n),  x2(n),  Xj(n)]T,  for  n  =1,2 N }  be 

the  J  correlated  image  scenes  which  contain  an  optical  signal  with  known  shape  and 
unknown  position.  Let  £  =  [  s(l),  s(2),  ....  s(N)  ]T  be  the  signal  pattern  N-  vector, 
and  &  =  [  bx  fo2,  •••■  by  ]T  be  a  /  -vector  of  signal  intensities  corresponding  to  the  / 
scenes  or  channels,  respectively.  The  two  hypotheses  which  the  adaptive  detector 
must  distinguish  are  given  by 

:  z(n)  =  z°(n) 

(I) 

Hx:  z{n)  =  £°(n)  +  &s(n) 


for  n  =  1,  2,  ,  N  where  <£°  is  the  vector  of  residual  clutter  noise-only  processes. 
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Thus,  under  hypothesis  H0  defined  in  Eq.  (1),  the  joint  probability  density  func¬ 
tion  of  the  Gaussian  clutter  vector  x(n )  is  given  by 


p(jc(n)  I  H0)  = 


1 


-~iT 


(2*) 


(j  +  Kyi 


I M  I 


1/2 


for  n  =  1,2 . N 


(2) 


where 


M  =  E  [(£<«)  -  E  £(n))(x(n)~  E  x(n))T  ]  (3) 

is  the  umcnown  covariance  matrix  of  z(n)  and  I M  I  *  0  is  its  determinant. 

It  was  demonstrated  previously  [3]  that  the  subtraction  of  a  space-  varying  local 
mean  from  the  image  can  yield  an  approximate  zero-mean,  near-white,  Gaussian  pr<  .- 
cess  with  a  slowly  space-varying  covariance  matrix.  Thus  for  N  sufficiently  small 
me  subimage  size,  N ,  can  be  chosen  so  that  matrix  M  is  approximately  a  constant. 
Experiments  it  dicate  that  such  residual  clutter  noise  is  also  approximately  indepen¬ 
dent  and  Gaussian  from  pixel  to  pixel  [3],  Therefore,  it  is  reasonable  to  assume  that 
the  residual  clutter  is  independent  from  spatial  sample  to  sample  (  see  appendix  A 
[3]). 

Next  define  J  xN  matrix  of  subimage  data  as  follows: 


X  =  U(1) . ,i(AO].  (4) 

The  j-th  row  of  matrix  X  is  x;( 2),  ...  Xj(N)]  for  j  =  1,  2,  J .  These  are 

the  N  observation  samples  of  the  y-th  image  scene  or  channel.  All  of  these  J  scenes 
or  channels  may  be  obtained  either  from  different  frequency  bands  of  the  same  image 
or  from  sequential  observations  of  the  same  scene.  Hence,  for  N  >  J  the  rows  of 
matrix  X  are  assumed  to  be  linearly  independent.  Then  because  of  the  mutual 
independence  of  the  components  of  X  in  Eq.  (4),  the  joint  probability  density  of  X 
under  H0  is  given  by 

N 

Po(X)  =p(x( l).i(2),  .  .  ,i(/V)  I  H0)=  FI  P&(n)  I  H0 ) 

n  - 1 
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(27r)N7/2IM  I*'2 


(5) 


In  terms  of  the  trace  function,  which  consists  of  the  sum  of  the  diagonal  entries 
of  a  square  matrix,  the  exponent  in  Eq.  (5)  can  be  re-expressed  as 


N  N 

S/(nWli(n)  =  Tr{M-]^x(n)rT{n))  . 

n= 1  n  =1 


Hence,  in  terms  of  (6),  p0(X)  in  (5)  becomes 

:±7><Af->£a;(«)ar(»)) 


(6) 


Po(X)  = 


(2nfm\M\Na 


(2n'f/J/2 1 M  ! Nf2 


(7) 


where 


M  =  XT  . 


(8) 


For  the  signal-plus-noise  hypothesis  Hx,  described  in  Eq.  (1),  one  computes  the 
conditional  mean  value  of  i(n),  given  hypothesis  H b  as  follows  : 


E[£(n)  I  //i]  =  EU°(*)  +  £  *(«)] 
=  b  s(n ) 


(9) 


for  n  -  \,  2,  ...  N .  Hence,  the  joint  probability  density  function  of  the  residual 
signal-plus-Gaussian  clutter  noise  vectors,  £(n)  for  n=l,  2,  .  .  .  N  under  hypothesis 
Hh  is  found  similarly  to  be 


Pi(X)  =  p(£( l),i(2),  &(N)  I  HO 


N 


-S 


=  FI  P(£(n)  !//))  = 


n  =1 


(2nfm\M  \Na 


(10) 
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where 


x  j,(n)  =  [x(n)  -  £  s(n)]  and  (11-1) 

Mll=-zZxll(n)zZ(n)  =  ±(.X-!£T)(X-,£T)r  .  (11-2) 

Nn= i  N 


Note  that  M  and  are  nonsingular  if  N  >  J,  (  e.g.  see  [5]  ). 

The  generalized  maximum  likelihood  ratio  test  (e.g.  see  [1])  is  shown  readily  to 
be 


|A£  {n/2  ^  k  ,  then  H , 

min..^  iN/2  <k,\henH0 
b  IA/*‘ 


Clearly,  the  tes:  in  Eq.  (12)  is  equivalent  to 


Z  = 


I M 


b 


>  c  ,  then  H  x 
<  c  ,  then  H0 


where  c  =  k2JN 


(12) 


(13) 


\M  1  =  I  XXT  I 


(14-1) 


and 


Itffcl  =  \(X  -b  S.T)(X  -  b  S.T)T\  ■ 


(14-2) 


Next  a  substitution  of  Eqs.  (14)  into  Eq.  (13)  yields 


_ 1  XX  T  I _ 

I  (X  -  b  S.T)(X  -  b  S.T)T  I 


>  c  ,  then  H  j 
<  c  ,  then  H0  ' 


or 


I 


:  xxT  i 


min  i  r 

b  1  Fk 


>  c  ,  then  H  j 
<  c  ,  then  H  0 


(15) 


where 
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F*  =  (X  -bST)(X  -b  ST)T  ■  (16) 

To  find  m^n  I  F^  I  one  can  start  by  expanding  Fb  as  follows: 

Fk  =  XXr  -  b  STXT  -  X(fe  ST)T  +  b  ST(h  ST)T 

=  XXT  -b  STXT  -XS  bT  +  b  bT(S.TS) .  (17) 


Then  by  adding  and  subtracting  the  term  simultaneously  in  Eq.  (17),  one 

can  decompose  the  J  x7  matrix  Fb  into  the  sum  of  three  matrices  in  such  a  manner 
that  only  one  of  them  contains  the  unknown  vector  b .  This  is  accomplished  by  first 
expressing  Eq.  (17)  as  follows: 


Fh  =  (StS) 


bb‘ 


bSTxT 

STS 


X$bT  (XS)(XS)T  (X$)(XS)T  ( XXT ) 
STS  (STS)2  (STS)2  STS 


(18) 


Hence,  by  a  "  completion  of  squares  "  operation  Eq.  (18)  becomes 


Fb  =  (STS) 


(£- 


-XLm-XL?  +  2°fl 
S.T  S.  sT  s  ST s 


(XS)(X$)T 
( STS )2 


(19) 


It  is  proved  in  Theorem  A.  2  of  Appendix  A  that  Fk  in  Eq.  (19)  is  positive 

VC 

=  Hence 
STS 

vc 

b  =  in  Eq.  (19),  one  obtains 
STS 


setting 


definite  matrix  for  any  vector  b  .  including  the  special  vector  b 


n  A  -  _  XX1  (XS)(XS)T 

8 TT  ars)2 


(20) 


where  B  is  positive  definite.  Thus  B  1  and  the  square  roots  B±V2  exist.  Note  that  B  is 
by  Eq.  (20)  independent  of  the  magnitude  of  vector  b  ■ 

To  calculate  the  determinant  i  Fb  I  one  can  establish  from  Eqs.(19)  and  (20) 
the  following  relationship: 


Fb  =(STS)B 


1/2 


+  B~m{b- 


--)(£~rrr)rF‘ 


1/2 


sTs  sTs 


B 


l.: 


(21) 
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Thus  the  the  deienninant  of  is 


Fk  I  =  (£T£/.IB  I  . 


lj  +  B-'u(h 


\n,u  \ru  M L\Td-1/2 


,  -)(&-  ,  . 
S  £  £  £ 


-)'S 


(22) 


Hence  finally,  the  denominator  of  Eq.  (15)  is  given  in  a  new  form  by 


min 


\Fk  \  =  (£r£)y  IS  I  **£ 


min 


h  +  S_l/2(& — —-)(£ — ^~)T  B~X!2 

£t£  £t£ 


(23) 


Since  B~xn(j2 — ^— )  is  a  Jxl-column  vector,  by  a  matrix  identity  proved  in 

£r£ 

Appendix  of  [3],  one  obtains 


lj  +  e~1/2(£— -££-)  (h-~-)TB'm 

£r£ 


£r£ 


-i+ia^L)r 


£r£ 


£r£ 


(24) 


where  B  1  is  positive  definite.  Therefore,  the  minimum  in  Eq.  (23)  is  obtained  when 
the  second  term  in  the  right  side  of  Eq.  (24)  vanishes.  But  this  happens  if  and  only 

if, 


Thus,  a  substitution  of  Eq.  (24)  into  Eq.  (23)  yields  finally, 
I Fkl  =  (£r£)y  IS  I 


=  (£r£/ 


\XT 

£r£ 


(X£)(X£  )T 

(£r£)2 


<YXr  (*£)(*£  )r  I 

(£'£t  I 


(25) 


(26) 


It  is  proved  in  Theorem  A.  2  of  Appendix  A  for  b  =  0  that  the  J  xj  matrix  XXT 
is  positive  definite,  so  that  the  square  roots  (XXTy'a  exist.  Hence,  the  denominator  of 
Eq.  (15)  can  be  expressed  further  as 


min 

b 


I Fb  I  =  IXXr  I 


h 


(XX T  y 1  l2(xs  )(xs  f  (xx T  r 1 a 

(£r£ ) 


(27) 
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Thus,  since  (XXr)~1/2(X£)  is  a  7xl  -column  vector,  Eq.  (27)  becomes 


aim  i  r-» 

t  IF‘ 


IXXrl (1 


(xs)r(xxTr'(xs) 

(S7^) 


(28) 


Finally,  a  substitution  of  Eq.  (28)  into  Eq.  (15)  yields 
j  >  r0  ,  then  H  [ 

1  _  (XS.)T(XXTy\Xl S)  <  ro  •  H o  {"9) 

LSr£) 

as  the  likelihood  ratio  test  function.  Clearly  by  Eq.  (15)  the  test  in  Eq.(  29)  is 
equivalent  to 


r 


(XS)T(XXTr'(XS)  "  r° 

(ir5)  <ro 


then  H  j 
then  H0 


(30) 


where  r  is  related  to  l  by  /  =  — - — ,  and  r0  =  1  -  — .  In  the  case  7  =  1  the  test  r 

1  —  r  c 

reduces  to  the  normalized  matched  filter  test, 


2ro-lhenWi 

(t’iXS’i)  < '« ■ mcn  H» 

a  well-known  CFAR  test  for  one  frequency  band  or  channel. 


III.  Detection  and  False  Alarm  Probabilities  of  Test 

In  order  to  find  the  probability  density  function  of  the  test  r  in  Eq.  (30)  on  both 
hypotheses  H 0  and  H  t.  one  partitions  matrix  X  as  follows: 

X  =  [  i(l)  I  i(2)  I  ...  I  iOV)]  (32) 

A 

where  z(n)  =  [j: jCzt )^c2(n )•  */0») I  is  the  7xl-column  vector  as  defined  in  Eq.  (1). 

By  (1)  and  (3)  one  has 

;  =  0,  I 


Cov  [  x(n  )  I  H,  1  =  M  for 


-9- 


Also 

E  [  i(n)  I  //0  ]  =  E  [  i°(n)  ]  =  0  and 
E  [  £(n )  I  H i  ]  =  E  [  £°(n )  +  bs (n)  ]  =  bs (n )  , 
or  in  terms  of  X  as  defined  in  Eq.  (32), 

E  [  X  I  H0  ]  =  0  and  E  [  X  I  H ,  !  =  b  ST  .  (33) 

Next  perform  a  whitening  procedure  on  £(n  )  by  defining 

z(n)  =  M~m  x(n)  ,  for  n  -  1,  2,  .  .  .  N  ,  (34-1) 

i  e  let 

Z  =  [i(  1)  I  1(2)  1  ...  I  z(N»  =  AT1'*  X  .  (34-2) 

The  whitening  procedure  in  Eq.  (34)  and  the  assumption,  that  the  residual  clutter 
samples  in  the  spatial  coordinates  are  mutually  independent,  produces  the  result, 

Cov  [  zt(m)Zj{n)  ]  =  5  (i-j  ,m-n)  (35) 

for  i,  j  =  \,  2,  ...  J  and  m,  n  =  1.  2,  .  .  .  N.  Here  5  (nsn  )  is  the  Kronecker  delta 


function  defined  by 

1  if  n  =  0  and  m  =  0 

5(n  jn  )  =  -  (36) 

.  0  otherwise  , 

and  z,(/n)  is  the  f-th  element  of  vector  z(m).  Then  by  Eqs.  (33)  to  (36), 

E  [  Z  l//0  )  =  0  ,  (37) 

E  [  Z  IW,  ]  =M~xa  b  £r  and  (38) 

Cov  [z(n)\Ht  I  =  l]  for  i  =  0,  1  .  (39) 


Evidently  by  the  transformation  in  the  Eq.  (34)  the  test  function  in  Eq.  (30) 
becomes  in  terms  of  Z  in  Eq.  (34-2)  the  expression, 

(ZS)T(ZZTr'(ZS)  '  r°  ’  thcn  H 1 
(£r£)  <  ro  >  fi100  H o 


r 


(40) 
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Now  let 

X  =  — 1  - 

(SrS)I/2  ' 

Then  the  test  function  in  Fq.  (40)  becomes,  using  Eqs  (40)  and  (41)  , 
r  =(Zi)r  (ZZrr‘(Zt) 


(41) 


(43) 


Also  by  Eq.  (41)  the  sum-of-squares  norm  of  i  is  given  by  llr.il  =  1  so  that  r  is  a  unit 
vector  in  the  "direction"  of  vector  £ . 


Now  consider  the  N  xN  orthonormal  matrix  U x,  which  carries  out  rotations  in 
N  -dimensional  space,  in  such  a  manner  that  unit  vector  i  is  transformed  into  the  new 
unit  vector, 

1  =  UXL  =  [1,0,  .  .  .  0]r  .  (44) 


Also,  let 

V  -  ZU\  -  [  v(l),  v(2),  .  .  y(/V)  J  ,  (45) 

Then  the  test  function  r  in  Eq.  (43)  reduces  to 

r  =  v(l)r(WrrIv(l)  .  (46) 

The  covariance  matrix  of  v(n),  for  n  =  1,  2,  ...  N  is  similar  to  that  of  z(n  ),  the 
only  change  of  the  statistics  of  the  v(«)  from  that  of  the  i(n)  is  their  mean  values 
under  hypothesis  Hx.  This  mean  is  derived  as  follows: 

E  [  V\HX  ]  =  E  [  V  UTl\Hl  ] 

=  M-'irz  t  LT  UTi  {S.T&)XU 

=  M~mt  (1,0,.  .  .  0 } (£t£)1/2 

=  [M~xa  b  (£r£)l/2  ,  Q - Q).  (47) 


From  Eqs.  (44)  and  (47)  a  figure  of  merit  or  what  might  be  termed,  the  generalized 
signal-to-noise  ratio  ( GSNR  )  of  the  test,  is  derived  as  follows: 
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( GSNR  )  =  E  [  yr(l)ltf ,  ]  E  [  v(l)IW,  ] 
=  (bT  M-lb) IIS  H2  =  a 


(48) 


Now  consider  a  further  simplification  of  the  test  function  r  in  Eq.  (46).  First 
separate  matrix  V  into  two  parts  in  such  a  manner  that 


N 


WT  =  V(1)VT(1)  +  £v(fi)vr(n)  =  v(l)vT(l)  +  Qt  , 

n~2 


(49-1) 


where 


s 


Q  =  Iv(rt)/(rt)  =  (v(2) . v(/V)][v(2),  ....  v(N)]t 

n-2 


(49-2) 


is  a  non-singular  J  xN  matrix. 

A  well-known  matrix  inversion  identity  applied  to  Eq.  (49-1)  produces  the 
result, 


(Wr)-‘  =  [v(l)vr(l)  +  Q  ]"'  =  [/  -  1  g~'  • 

1  +yr(l)Q-1v(l) 


(30) 


A  substitution  of  Eq.  (50)  into  the  test  function  r  in  Eq.  (46)  yields  the  test  functio 
r  as  the  new  expression. 


r  _  ^(DQ-'vd)  ri 

1  +  vr(l)2_,v(l)  1+r, 


(51) 


where 


rx  =  yr(lX2"1v(l) 


(52) 


It  is  desired  now  to  find  the  probability  density,  /(r,!//,)  ,  of  rx  in  Eq.  (  52). 
First  by  Eq.  (49-2),  define  D  =  [v(2),  .  .  .  ,  v(/V))  and  re-express  Eq.  (52)  in  the  form. 


■  2/  V  (1) 


rr  i  y(i) 


r ,  =  II v (1) II (DD1  ) 

1  -  ii  v  ( i )  it  ii  v  ( n  i 


) 


(54) 


Then  normalize  the  J  -component  vector  v(l)  as  follows  : 
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t  _  Yd) 

**  II  v  (1)11 

Hence  by  Eq.  (55)  one  obtains  r  in  Eq.  (54)  in  the  form, 
r i  —  Hv(l)ll2(^r(DDrr1^)  =  Hv(l)ll2e  , 

where 

e  =  (^(DDVD- 


(55) 

(56-1) 

(56-2) 


Now  one  can  farther  process  the  term  e  in  Eq.  (56)  by  conditioning  on  the  ele¬ 
ments  of  v(l)  so  that  £  can  be  treated  as  a  normalized  constant  vector.  Then  since  £ 
has  unity  magnitude,  there  exists  a  J  x7  orthonormal  matrix  U2  such  that, 

C/2£=[1,0,  ...0]r.  (57) 

Next  apply  this  transformation  to  matrix  D ,  defined  before  Eq.  (54),  by  letting 

H=U2D  =U2[v(2),  ...,v(N)].  (58) 

Then  the  term  e  in  Eq.  (56)  has  the  simple  form, 

e  =  ZT(DDTr'i={  1,0 _ OKtfZ/Vn.O _ 0]r  .  (59) 

Clearly  H  in  Eq.  (58)  has  the  exactly  the  same  statistical  properties  as  D ,  under  the 
assumption  that  v(l)  is  given. 

Now  partition  H  as  follows: 

hi 

H  =  ,  (60-1) 

Hb 

where  hA  is  the  /V-l-column  vector  and  HB  is  the  (Z-l)xJV-l)  matrix.  Then  by  Eq. 
(59), 

hi  hA  hIHI  *  raa  Rab 
o hht )-‘  =  = 

HbIIa  hb  hb  [Rba  rbb 


(60-2) 
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According  to  the  Frobenius  relations,  e.g.  see  [8]  or  [9],  for  a  partitioned  matrix, 

rm  = [bib  -  r>  =  [tin  -  htb(hbhtbt'hb) ^  r1 

1  1 


hid hlPxhA 
A  substitution  of  Eqs.  (60)  and  (61)  into  Eq.  (59)  yields 
1 


(61) 


e  - 


hi?  itlA 


(62) 


where  Px  =  IN_i  -  HB(HBHByxHB  is  a  projection  operator  such  that  P,2  =  Px  and 
Tr  (Pj)  =  N-J .  In  the  same  manner  used  for  the  projection  matrix  P  in  [1]  it  is  not 
difficult  to  show  that  Pt  has  N-J  unity  eigenvalues  and  J-\  zero  eigenvalues.  Thus 
P  x  can  be  diagonalized  to  the  form  , 


Ul  P,  t/3  =  A,= 


l N-J 

0 


0 


(63) 


where  IN~j  is  the  (N-J)x(N-J)  identity  matrix.  By  arguments  similar  to  those  used 
previously  in  [1],  one  finds  also,  under  the  assumption  that  v(l)  and  Px  are  given, 
that 


where 


N-J 

hTAP\hA  =nr  n=  Eh,2. 

,=i 


n  =  hlUj  a1/2 


(64-1) 


(64-2) 


is  a  1  x(/V-/)-column  vector.  The  conditional  joint  probability  density  function  of  n 
is  subject  to  the  normal  density  function,  N  (  Q  ,  lN_j  ),  i.e. 


Prx (hi-  ■  ■  ■  .  h(N-y>  I  v(l).f> \)  -  N (QJn-j) 


(65) 


Since  pn(hi,  •  ■  •  .  h(/v-/)  I  v(l).^i)  in  Eq.  (65)  does  not  depend  on  v(l)  and  Pb 
the  vector  y(l)  and  matrix  Px  must  be  statistically  independent  of  vector  r\.  Also  by 
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Eq.  (62)  one  obtains  the  ratio  r,  in  the  form, 

r  1) 

nn  Xn? 

i-I 


(66) 


in  term  of  the  magnitudes  of  these  two  vectors  only.  Thus  by  the  independence  of 
vectors  n  and  v(l)  and  a  use  of  the  Corollary  2  in  [7,  pp  52],  one  has  the  probability 
density  function. 


/(r,lH,)  = 


J-  2  — 

T~  2 
r  i  e 


S(^’i/)(I  +  ri)2 


N  \  or  j 
—  ifiC- -.7/ .  ,). 


2  2  2(1  +  r,) 


(67) 


of  r,  in  Eq.  (51)  under  hypothesis 

Finally  by  using  the  relationship  of  r,  to  r  in  Eq.  (51)  the  probability  density 
function  of  the  test  function  r  under  hypothesis  Hx  is  given  by 


f(r\H  d  = 


rtf) 


N-J-  2  3-2 

2  ,  2 


rc 


(1  -r)  2  r  2  e  2  iF{(— ) 

N-J^J,  1  1  2  ’2  2 


(68) 


2  ■)F(t) 


for  0  <  r  ^1  where  a  is  the  generalized  SNR  in  Eq.  (48)  and  F}(a;b  ;x)  is  the 
confluent  hypergeomitric  function.  Clearly,  if  no  signal  is  present,  then  a  -  0  .  Thus 
Eq.  (67)  reduces  in  the  H0  hypothesis  to  a  Beta-function  probability  density  of  form. 


f(r\H0)  = 


rff) 


N-J  J 

R-^)r(y) 


N-J -2  J-  2 
(1  -r)  2  r  2 


for  r  <,  1 


(69) 


Finally  in  terms  of  the  above  probability  density  functions  in  Eqs.  (69)  and  (68) 
the  probability  of  a  false  alarm  is  found  by 
1 


Pfa  =  ]f{r\H0)dr 


(70) 


-15- 


and  the  probability  of  detection  by 

<30 

PD  =  jf(r\H  x)dr  , 
ra 


(71) 


IV.  Performance  Analysis 

Performance  curves  of  the  probability  of  detection  in  Eq.  (71)  versus  SNR  for  a 
given  false  alarm  probability  with  respect  to  various  values  of  N  and  7  are  computed 
in  this  section.  First,  in  Fig.  1  the  probability  of  detection  is  calculated  as  a  function 
of  the  generalized  signal- to-noise  ratio  (GSNR)  a  in  Eq.  (48)  for  several  different 
values  of  parameter  N.  These  curves  demonstrate  that,  for  a  fixed  GSNR,  a,  the 
CFAR  detector  has  a  higher  detection  probability  if  more  samples  are  used. 

In  order  to  compare  the  detection  performance  improvement  for  different 
numbers  of  signal -plus-noise  bands,  the  probability  of  detection  in  the  single  scene  of 
maximum  SNR  is  compared  with  the  probability  of  detection  in  two  correlated 
scenes,  i.e.  for  7=2.  To  accomplish  this,  the  GSNR  in  7  correlated  scenes  is  related 
to  the  maximum  SNR  of  the  7  scenes.  Consider  first  the  case  of  7  =  2. 

For  7=2  let  the  maximum  signal-to-noise  ratio  in  the  2  correlated  scenes  be 
given  by 

.  b ,2  IIS  II2 


a  = 

<*\ 

Then  the  GSNR  in  Eq.  (48)  for  7  =2  is 
a  =[  bx  b2]M~l[b{  b2f  \\  S»2  , 

where 


(72) 


(73) 


M  = 


of 

o2i 


a,22 

of 


(74) 


-16- 


Hence, 


A/-1  = 


0,2022  -  Ojj  [  "°21 
and  as  a  consequence  one  has 


-C{2 

o? 


(75) 


t  I,  Cl, 2  (62<32  +  6 2 a2  -  26 1  62  0,2)  IIS 


a  =  [ b i  biY  « £  »*  = 


_2_2  _2 
0i02  -  0)2 


(76) 


1  6  2  II  £  II  6£  115  II  26 1^2cr12 


1  - 


_2 
°12 
2  2 
<*f<*2 


■( 


Of 


02 


2^2 


ora 


1  u2 


a 


1  - 


62  °i2  62  Oi 

— 0  +  “I "7  “  2-T—L>- 

a,2  bf  a{  6,  o2 


(77) 


a2a2 


In  terms  of  the  normalized  correlation  coefficient, 

0,2 


0,02 


(78) 


and  the  ratio, 

X2^ 


2  A  (smh  =  biwMol 

(SNR) ,  6  2  II 5  II 2/o2  ’ 


(79) 


of  SNR  in  the  second  scene  to  SNR  in  the  primary  scene  of  maximum  SNR,  Eq.  (77) 
can  be  re-expressed  as  follows: 

a=Ga'  (80- n 

where 

G  =  (1  +  ■-  (80-2) 

1  -p2 

is  the  gain  in  signal-to-noise  ratio  of  a  detector  which  uses  7  =  2  scenes  over  a  single 
detector  which  uses  that  scene  with  the  maximum  SNR.  The  following  lemma  is  now 
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proved. 

Lemma  1: 


q  _  (1  +  X2  -  2 Xp)  >  ! 

l  -  p2 


for  all  X  and  p  with  0  <  X,  p  <  1  such  that  condition  X  =  1  and  p  =  1  is  not  true. 


Proof :  Start  with  inequality  (  X  -  p  )2  >  0,  with  equality  if  and  only  if  X  -  p.  This 
implies  X2  -  2Xp  >  -p2  which  in  term  yields  1  +  X2  -  2Xp  £  1  -  p2.  Thus  if  X  *  p, 

then  by  above  inequality,  G  =  +  >  1  and  Lemma  is  true.  However  if 

1  -  p2 


X  =  p  and  0  <p  <  1,  G  =  — 2Sil  -  ]  and  again  Lemma  is  true.  Finally  if 

X  =  p  and  p  =  1.  Clearly  for  this  final  case  G  is  indeterminate.  Thus  Lemma  1  is 
proved.  The  above  Lemma  is  generalized  to  any  positive  integer  /  in  Appendix  B. 


The  above  Lemma  shows  that  the  new  CFAR  detector  algorithm  of  a  target  in 
two  correlated  scenes  is  always  better  than  a  CFAR  detector  in  the  scene  with  max¬ 
imum  SNR.  Fig.  2  illustrates  the  probability  of  detection  for  a  false  alarm  probability 
of  PFA  =  1(T5  as  a  function  of  a'  for  J  =  1.  2,  N  =  49,  X  =  0.2  and  p  =  0.95,  i.e.  for 

G  =  6.77.  This  shows  that  this  detector  using  2  scenes  with  ( SNR)2  =  ~(SNR)\  has  an 

approximate  8.5  dB  SNR  improvement  over  a  detector  which  uses  the  single  scene 
with  maximum  SNR. 


A  comparison  of  the  SNR  in  Eq.  (80-1)  is  now  made  with  Eq.  (46)  in  [1]  for 
K  =  1  the  single  noise-only  reference  scene  case  in  [1].  By  Eq.  (46)  in  [1]  the  SNR 
for  a  scene  with  signal,  using  one  reference  without  signal,  is  given  by 

a  =  (Cy/Oy\x)a'  =  — — rr  =  Ga’  (81) 

(1  -  P  J 


in  terms  of  p  and  a',  the  SNR  in  the  scene  with  signal.  By  Eq.  (80-1)  the  GSNR 
depends  on  p,  the  correlation  coefficient,  but  also  on  X,  the  ratio  of  SNR’s  in  both 
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scenes.  It  is  evident  from  Eqs.  (80-1)  and  (81),  that  as  long  as  the  inequality 
1 

—  <  p  <  1  holds,  the  gain  in  SNR  in  Eq.  (80-1)  is  always  less  than  or  equal  to  the 
gain  in  the  SNR  given  in  Eq.  (81). 

Using  the  correlation  coefficient  p  =  0.81,  a  computational  comparison  is  made 
in  Fig.  3  with  the  results  given  in  11]  for  K  =  1.  These  curves  in  Fig.  3  illustrate  that 
for  the  same  p  the  probability  of  detection  curves  for  the  new  CFAR  detector  are 
limited  to  the  probability  of  detection  curves  given  in  [1]  for  K  =  1.  The  leftmost 
curve  in  Fig.  3  is  as  derived  under  the  assumption  that  there  is  no  signal  in  the  refer¬ 
ence  scene. 

Figs.  4(a)  and  (b)  show  typical  32  x  32  subimages  in  two  different  optical  bands 
(  the  green  and  red  bands  )  of  the  San  Diego  area.  A  5x5  signal  with  the  pattern, 
given  in  Fig.  5,  is  implanted  in  both  of  these  green  and  red  images  with 
(SNR  )t  =  0.37  and  (SNR  )2  =  0.2x0.37.  A  local  mean  is  subtracted  from  both  of  these 
subimages.  The  resulting  residual  images  are  approximately  zero-mean  Gaussian 
processes.  The  CFAR  test  given  in  Eq.  (30)  is  calculated  for  each  pixel.  The  test 
statistic  of  this  CFAR  test  is  plotted  pixel  by  pixel  in  Fig.  4(c).  A  target  is  detected 
with  a  threshold  determined  by  a  Pfa  =  10~5  .  The  results  of  this  test  are  illustrated  in 
Fig.  4(d)  which  shows  that  the  target  was,  in  fact,  detected. 

0  0  1  1  f 
0  0  110 
0  0  10  0 
0  110  0 
10  10  0 


Fig.  5  Target  signal  template 
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In  order  to  demonstrate  the  theoretical  SNR  improvement  in  Eq.  (80-1),  a  com¬ 
puter  simulation  was  performed  to  determine  the  required  SNR  in  the  primary  scene 
of  maximum  SNR  needed  to  detect  a  target  in  single  primary  scene  and  in  J  -  2 
scenes.  A  computer  simulation  was  made  similar  to  that  developed  in  [1].  The  results 
of  this  new  CFAR  test  are  shown  in  Table  1,  where  k  =  0.2.  Using  5  different  subim¬ 
ages  the  average  improvement  of  a  detector  using  J  =  2  scenes  over  the  single  pri¬ 
mary  scene  is  6.21  dB. 


V.  Conclusions: 

Under  the  same  assumptions  for  optical  noise  clutter  used  in  [1],  a  generalized 
constant-false-alarm-rate  (CFAR)  algorithm  is  developed  for  detecting  the  presence  of 
an  optical  signal  of  non-zero  intensity  in  J  signal-plus-noise  bands  or  channels.  For 
many  applications  this  new  algorithm  is  more  flexible  and  practical  than  the  one 
given  in  [1],  If  /  =  1  ,  the  resulting  test  reduces  to  the  standard  normalized  matched 
filter  test  for  finding  a  signal  in  clutter  of  unknown  and  varying  intensity. 

Both  theoretical  and  computer  simulation  results  show  that  the  SNR  improve¬ 
ment  gain  of  this  new  algorithm  using  multiple  band  scenes  over  the  single  scene  of 
maximum  SNR  is  always  greater  than  one  and  in  some  .ases  it  can  be  substantial.  A 
comparison  of  SNR  gain  between  this  new  detection  algorithm  and  the  one  given  in 
[1]  illustrates,  that  for  the  same  correlation  coefficient  p  of  related  scenes  the  proba¬ 
bility  of  detection  curves  for  the  new  CFAR  detector  are  limited  to  the  probability  of 
detection  curves  given  in  [1]. 
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CFAR  Algorithm  Simulation  Results  for  the  Optical  Color  Image  with 
N  =  49  Samples  for  Covariance  Estimation 


Subimage 

Location 

J=  1 
(dB) 

J  =  2 
(dB) 

(100,  10) 

6.67 

0.96 

( 234,  30 ) 

9.34 

2.04 

(  50,  310) 

7.32 

5.18 

( 200,  200 ) 

2.30 

-4.32 

(65,  1) 

3.80 

-4.20 

Average 

5.88 

-0.33 

Improvement 

Factor 

— 

6.21 

Table  1 
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Appendix  A 


Theorem  A.  i: 

A  is  nonnegative  definite  if  and  only  if  there  exist  a  matrix  W,  such  that 
A  =  WWT  .  For  a  proof,  see  [5,  pp  257]. 


Theorem  A.  2: 

The  J  xj  matrix  Fb  =  (X  -  b  S.)(X  -  b  S.)T  in  Eq.  (16)  is  positive  definite  for 
any  arbitrary  vector  b  . 


Proof  :  Since 

Fb  =  (X  -b  Sr)(X  -  b  ST)T  (A-2) 

=  WWT 

where  W  =  (X  -  b  &r)-  By  using  Theorem  A.  I  matrix  Fk  is  nonnegative  definite, 
and  all  the  eigenvalues  of  satisfy  X}  >  0  .for  ;  =  1,  2,  ...  J.  But  by  previously 
showing  that  were  nonsingular,  F^  is  also  nonsingular  by  Eq.  (14-2).  This  means 
that 


det  (Fk)  =  .  .  .  Xj  *  0  . 


or  Xj  >  0  for =  1,  2,  .  .  .  J.  Hence  Fb  is  a  positive  definite  matrix.  In  particular  if 
XS 

one  chooses  b  =  then  by  Eq.  (26) 


B 


XX T  _  (XS  )(XS  )T 
S.TS.  (£rS)2 


(A-3) 


is  a  positive  definite  matrix. 
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Appendix  B 

In  this  Appendix,  Lemma  1  of  Sec.  IV  is  generalized  to  an  arbitrary  positive 
number  J.  By  Eq.  (48)  the  generalized  signal-to-noise  ratio  (GSNR)  for  J  bands  is 
given  by 

a=[blb2 . bJ]\ri[blb2 _ bj  ]r  II £  H 2  (B.l) 

where 


M  =  Ei 


' 

r 

r  i 

r 'i 

< 

£ 

\ 

H 

£ 

i 

H 

1 

> _ 

*1*1 
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and 


A  =  [x2  x3,  .  .  .  ,  Xj  ]r  . 


In  terms  of  the  above  notation  the  inverse  of  M  in  Eq.  (B.2)  can  be  put  in  the  form, 

(B.3) 
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By  the  Frobenius  relations  in  [8]  or  [9],  one  has  the  results 
A  =  {K.  x  -Kx  xK~lK„  r1  = - ! - r -  , 

JT  jX  j  ii  Sll/  r*  rs  r/  —J  w- 

B  =  ~(KX[X[  -  KXiXK^Ksx)-'KXiXK;1  = 
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D  =  Kxl  +  , 
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(B.4-1) 


(B.4-2) 


(B.4-3) 


(B.4-4) 


for  the  submatrices  defined  in  Eq.  (B.3).  Hence,  one  can  rewrite  Eq.  (B.3)  as 


A/-1  =  — 


~KxlKIlX 


K~lA  +  K-lK&Kx^K-\ 


(B.5) 
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where 

A  =  (**,*,  -  (B.6) 

Note  that  A  >  0  since  M  is  the  positive  definite  matrix  (  see  Sec.  I  Eq.  (3).  Next  by  a 
substitution  of  Eq.  (B.5)  into  Eq.  (B.l),  one  obtains 

a  =  [bxb2 - bj]M~'[bxb2 - bj  ]r  II  £  li 2 


=  [b\  bTW~\bx  bT  l7  US  H2  (B.7) 

=  Mb ?  -  bxhTK-lK^x  -  KXllK~lbb ,  +  bTK~lb  A  +  t>T K^K^K^b  }  IIS  II2  , 


where 

b.  =  [  *2,*3 . bj]T  . 


In  Eq.  (B.7)  define  the  general  normalized  correlate  coefficient 
V-  ~  *SS  A5Ct 


(B.8) 


and  the  ratio 


A  bTK~lb 
A  4  — - ; — T- 


(B.9) 


of  the  GSNR  in  the  other  scenes  to  SNR  in  the  primary  scene  of  maximum  SNR. 
Then  Eq.  (B.7)  becomes 

a=Ga’  (B.10) 

where  G  is  the  generalized  gain,  given  by 


G  1  -  2lTo.  ±  -  ar&) +  &ra)2 

i-pre 


(B.l  1) 


Note  by  Eqs.  (B.6)  and  (B.8)  and  the  comment  following  Eq.  (B.6)  that  prp  >  1. 
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Lemma  :  Let  G  be  the  generalized  gain  function  defined  in  Eq.  (B.10),  then  G>1. 
proof  :  Note  first  that 

G  >  1  iff  1  -  2&rg  +  ^r^(l  -  pTp)  +  Qj Q)2  £  1  -  prp. 

iff  eTa  -  2 hTo.  +  -  ero)  +  &rco2  ^  o  (B.12) 

where  "  iff  "  denoted  "  if  and  only  if  ",  the  logical  equivalence.  Next  let 

Qr0_  =  a2  (B.13-1) 

&r2,  =  P2  (B.  1 3-2) 

&rp  =  apcosy  (B.  1 3-3) 

then  using  (B.13-1)  to  (B.13-3),  the  last  statement  in  (B.12)  becomes 

G  >  1  iff  a2  -  2aPcosy  +  p2(l  -  a2)  +  a^cos^  >  0  .  (B.14) 

There  are  two  cases  to  consider, 
case  I  :  If  p2  >  1,  then  since  1  -  prp.  >  0,  one  has 

a2  -  2aPcosy  +  p2(l  -  a2)  +  a2p2cos2y  S  a2  -  2aPcosy  +  (1  -  a2)  +  a^cos^  (B.  15) 

=  1  -  2apcosy  +  a2p2cos2y  =  (1  -  apcosy)2  >  0  . 

case  II  :  If  p2  <  1 ,  then 

a2  -  2  apcosy  +  p2(l  -  a2)  +  a2p2cos2y  =  a2  -  2aPcosy  +  p2  -  a2p2(l  -  ccs2y) 

>  a2  -  2aPcosy  +  p2  -  a2(l  -  cos2y)  (B.16) 

=  (acosy  -  P)2  >  0  . 

Thus  G  >  1  and  Lemma  is  proved. 
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